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Abstract 



We construct a semiclassical expression for the Husimi function of au- 
tonomous systems in one degree of freedom, by smoothing with a Gaussian 
function an expression that captures the essential features of the Wigner func- 
tion in the semiclassical limit. Our approximation reveals the "center and 
chord^ structure that the Husimi function inherits from the Wigner function, 
down to very shallow valleys, where lie the Husimi zeroes. This explanation 
for the distribution of zeroes along curves relies on the geometry of the classi- 
cal torus, rather than the complex analytic properties of the WKB method in 
the Bargmann representation. We evaluate the zeroes for several examples. 
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I. INTRODUCTION 



The features that distinguish integrable from chaotic motion in classical mechanics mani- 
fest themselves most clearly in phase space. This is one of the reasons for the great interest in 
the so called quasiprobability distribution functions" in phase space within the semiclassical 
theory of quantum states. These distributions are defined as the symbols associated with 
the density operator p in some representation of quantum operators PJ^. It is expected 
that these representations of quantum states show the differences between an integrable, 
or a chaotic classical counterpart in the semiclassical limit {h ^ 0). Among these phase 
space representations of quantum states, the Wigner function {i.e. the symbol of the den- 
sity operator p in the Weyl representation) and its smoothing by a Gaussian function, the 
Husimi function, are of paramount importance. In fact. Berry showed that the peak of 
the amplitude of the Wigner function is located very close to the curve of constant energy, 
for a pure state of an autonomous system with one degree of freedom, collapsing onto a zero- 
width distribution {i.e., a delta function) over that curve in the classical limit {h = 0) (see 
also [@-|^). Ozorio de Almeida and Hannay generalized this picture for states supported by 
invariant tori of higher dimensions . In all such systems, the semiclassical analysis of the 
Wigner function reveals an interesting geometrical structure of chords and centers" 
that determines the phase of the oscillations of the Wigner function as the point x = {q, p) is 
varied within the torus. This phase is proportional to the symplectic area -or center action- 
bounded by the torus and the chord, centered on x, joining two points of the torus. The links 
of this ^^semiclassical geometry" to the generating function formalism of classical mechanics 
and the path integrals of quantum mechanics are reviewed in reference 

Although the oscillations of the Wigner function thus reflect legitimate structures of 
classical mechanics, its positive definite Gaussian smoothing, the Husimi function, is much 
closer to a classical Liouville density. The density peaks near the region of classical motion, 
decaying exponentially in classically forbidden regions |]14|,|15|. The first impression is that 



smoothing cancels all trace of the centre and chord skeleton of the Wigner function. However, 
we shall show that very delicate effects are still discernible. 

First, we must recall that the Husimi function can also be viewed as the mean value 
of the density operator p in coherent states, whose holomorphic (entire) part is called the 
Bargmann function 0. This function corresponds to the wave function for the quantum 
state in a representation of the quantum mechanics introduced by Bargmann (in the 



case of the Heisenberg-Weyl group), where the basis for the Hilbert space is made of coherent 
states not normalized to unity. Thus, in the Bargmann representation the wave functions 
are holomorphic (entire) functions of the variable z = -^{(3q — ip/ (3), acting as a phase space 
coordinate. The analiticity of the Bargmann function compels its zeroes and those of the 
Husimi function to be isolated for 1-D systems. Leboeuf and Voros |16| have shown that in 
many cases the distribution of Husimi zeroes is completely different for chaotic maps, where 
they are spread out, as opposed to integrable maps, where they are distributed along curves. 
Only the latter alternative is available for systems with continous time and one degree of 
freedom. Furthermore, these lines of zeroes cannot occur close to the energy shell where 
the smoothed Wigner function has a non-oscillatory peak. The lines supporting zeroes may 
only be found in regions where the Husimi function is already exponentially small. 

In these circunstances, we can only expect to predict the general pattern of zeroes with 
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a very delicate ^^subdominanf semiclassical theory. This is the case of WKB type of theory 
developed by Voros for the Bargmann representation [^J, which predicts zeroes on the anti- 
stokes lines where two or more branches of the complex action have the same amplitude. 
The zeroes along these lines are selected by the condition that the imaginary part of the 
complex action be an integer multiple of vr. However, this approach has practical difficulties 
to obtain explicit formulae even to leading order in h. First, it is generally very difficult 
to obtain analitically the branches of the classical energy curve in complex coordinates, in 
order to calculate explicitily the complex action [i.e. the phase in the WKB wave function). 
Second, an approximation valid anywhere outside the neighborhood of the energy curve, 
requires the analytic continuation of the functions that define the branches; this needs the 
analyticity of the Weyl symbol, Hw, for the quantum Hamiltonian. 

Let us thus return to the picture of the Husimi function as a smoothing of the Wigner 
function. Since our approximation does not maintain explicitily the analytical properties, 
one could not expect to establish that there exist isolated zeroes in this way, but we can seek 
for shallow valleys, even in the region where the Husimi function is already exponentially 
small, and for oscillations along their bottom as indications of where the zeroes may lie. The 
simplest guess is that the two dominant regions in the evaluation of the Husimi function are 
the neighborhood of the centre of the Gaussian and the maximum of the Wigner function 
along the energy curve S. Since we know that a zero will only be found when the Gaussian is 
far removed from the energy curve, we use Berry's simple cosine-oscillatory representation of 
the Wigner function for the local approximation. For the contribution of the region near the 
energy curve, we start from the even cruder classical aproximation that the Wigner function 
is a ^-function along S. After the smoothing, the first term remains cosine-oscillatory with 
essentially the same phase (proportional to the center action, upon small corrections), but 
now damped by an exponential function which decreases, essentially with the length of the 
chords (upon small corrections). The second term is everywhere positive, smooth and peaked 
in the energy curve S. The combination of both produces a positive smooth expression, 
peaked along the curve S oscillating in a valley of local minima that approach the zeroes of 
the Husimi function when ^ — »• 0. This expression is valid only inside the energy curve and 
depends only on the properties of the torus S. 

The paper is organized as follows. In Section |I| we summarize important results con- 
cerning Wigner and Husimi functions. In Section ^TT| we find a semiclassical expression of 
the Husimi function for the case of a particle in a box, as a simple model of the geometrical 
approach. In Section |rV| we introduce our geometrical approach to the distribution of the 
Husimi zeroes in 1-D systems. In Section we apply this approach to the problem of a 
particle under the action of a constant force. This example corresponds to an unbounded 
problem whose convex energy curve is open. Finally, in Section ^ we present the results 
for the case of a particle subject to an asymmetric anharmonic potencial, as example of a 
general system with a convex and closed energy curve. 

II. REVIEW OF WIGNER AND HUSIMI FUNCTIONS 

Quasiprobability distribution functions are symbols associated with the density operator, 
p, in some representation of quantum operators The Wigner function is the Weyl 
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symbol of the density operator. The symbol of an operator A, in the Weyl representation is 
given by the function, 

Aw{^) = Jd^g (g + e,/2|A|g-e,/2) exp H^pA] , (1) 

where x = {q,p) (all integrations in this work run from — oo to +00 unless indicated). So, 
in the case of pure states in 1-D systems, the Wigner function is 

^^""^ ^ (i) P"^^""^ = (i) / "^^^ + ^./2|^)(^k - e./2) exp HpUn] . (2) 

We remark that since Tr[p] = / ^pvk(x) = 1 for normalizable states \ip), while diverging 
otherwise, the prefactor in (0) will not be considered for unbounded states. 

The semiclassical analysis of this function was first developed by Berry for the case of 
an eigenstate of energy E, in nonrelativistic 1-D systems, where the classical Hamiltonian 
is of the form 

H(x) = pV2m + V(g) (3) 

and the "torus" is the smooth convex curve, S, of constant energy (H(x) = E). We briefly 
summarize the results in which are important for this work (for more details see also 

ID- 

a) The simple semiclassical approximation. 

This is obtained by replacing the primitive WKB functions {i.e. the semiclassical solution 
of the time independent Schrodinger equation) in (0) and evaluating the integral by the 
stationary phase method. The result is symmetric in q and p and depends only on the 
geometry of the classical curve S, 

Wscl{^) = V cos 1^ - ^1 , (4) 

where the function S'j(x) is the symplectic area bounded by the energy curve S and the 
chord C,, centered in x, that joins the points x+ and x_ on the "torus" (see FIG.[l|). The sum 
is over all the chords centered on x and uj is the frequency of the classical motion around S. 
The Dj{x) are the skew products of the phase space velocities at the tips of the chord., 

Dj(x) =x_ Ax+=p_g+-g_]3+ , (5) 

representing the area of the parallelogram formed by the pair of vectors (the point x_|_ is 
reached after the point x_ in the classical motion along £, see FIClH). Outside of the convex 
energy curve £ there are no chords, so Wschi?^) = 0. 

The Wigner caustic labeled as C in FIG.|I], is the border of regions with different numbers 
of chords: within the caustic there are three chords, on the caustic two and outside only one. 
On the Wigner caustic and on generically two pairs of stationary points coalesce; in both 
cases the simple method of stationary phase is inapplicable. We stress that two of the three 
terms in (^) diverge as a smooth side of C is approached from the inside (FIG.|I|). The 



4 



reason is that the phase space velocities, x+ and x_, are parallel, so the area Dj{x.) is zero. 
We also remark that when the eigenstate is normalizable, the prefactor in (^, that arises 
from the correct normalization of the primitive WKB function, does not give the correct 
normalization of the Wigner function, 
b) The uniform approximation. 

Simultaneous consideration of a pairs of stationary points in (|^) yields 



Wscl{^) 



V2 



E 



35, (x) 



1/6 



Ai 



3gj(x) 
2h 



1 2/3' 



(6) 



This is an uniformily valid approximation not only as x moves onto S, but also when x lies 
on the convex side of S (H(x) > E) where the stationary values and the function ^(x) are 
imaginary. However, further refinements of the method of stationary phase are required to 
obtain an approximation uniformly valid over C. 

On the concave side of £ (H(x) < E) and not too close to S, Six) is large in comparison 
with h so the Airy function can be replaced by its asymptotic form for large negative 
argument [0, 



1 f 2 

Ai{— w} ^ —=[w]~^^'^cos\ -[w 
y/n 13 
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(7) 



and (^) is recovered. On the convex side of S the Airy function has positive argument and 
the semiclassical Wigner function decays exponentially away from S. 

When the eigenstate is normalizable, the uniform Wigner function is correctly nor- 
malized. 

c) The transitional approximation. 

Very close to £ an expansion of (0) yields. 



Wscl{^) 



1 



where 



271 u 



-1} 



Ai 



[rS(x)]l/3 1^ [^^S(x)]l/3 



[H(x) - E] 



Hpp + Hp Hgg + 2 H500 H„ H 



q ' 



(8) 



(9) 



with all the partial derivatives of H evaluated at x. B(x.) remains finite as x moves onto S. 
d) The classical limit. 

This corresponds to the limit when ^ = and is obtained by letting h ^ in (§) and 
using the result. 



lim - Ai I — 



5(y) 



(10) 



to give. 



— <5[H(x) - E] . 



(11) 
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Along the Wigner caustic the modulus of the Wigner function takes large values. However, 
the infinitely rapid oscillations along C cancel the amplitude of the delta function. In the 
case of normalizable eigenstates, (^) and (|1T]) are correctly normalized. 

The prefactor in @ and the normalization constants in (^), and (pJ]) , are defined 
for normalizable states. When the states are not normalizable, the formulae are still valid, 
but now it is possible to define the normalization using the orthogonality conditions if the 
wave function belongs to an orthogonal set 0J§]. For these cases, the classical frequency a; 
should not be included in the formulae. An example of this type of normalization appears 
in Section |V|. 

From the fundamental quasiprobability" property (see e.g. 0), 

{^\A\^) = jd^ Aw{^) IV^(x) , (12) 
we obtain that the "scalar product" of Wigner functions, 

j dxiy^(x) W^{^) = {J-^ 1(01^) |2 , (13) 
is always positive definite, including the projection onto positions, 

j dpW^{^) = \{qm\ (14) 



giving the probability density (see pO| ). 

Recently, Ozorio de Almeida 0] reviewed the link between the semiclassical geometrt/^ 
underlying the Wigner function and the generating function formalism of classical mechanics. 
This is based on canonical conjugate variables, the centers, x = (g,p), and the chords, 
^ = (Cq,Cp), as an alternative in the description of the classical evolution. Hence, instead of 
specifying the 2 L- dimensional initial, x_, and final, x+, points in phase space, we can give 
the vector ^, joining each pair, and the position x of its center. The evolution is described 
by the canonical transformations given by the center or chord generating functions [i.e. the 
center or chord actions respectively), from which we obtain the corresponding canonical 
conjugate variables by different at ion: 

«x)^-^ , Ux,^^. (15) 

The center action S'(x), for fixed energy, is the function in the expressions (^) and (^. 

The Husimi function is another quasiprobability distribution function; it corresponds to 
the normal symbol of the density operator in the diagonal coherent states representation |Q] 
(or Husimi representation [|l^]). In this representation the normal symbol of an operator A 
is the expectation value, 

A^(X) = {Q^\A\Q^) , (16) 



where the \^x) are the minimum uncertainty states known as coherent states [|TT],|12[. These 
states are eigenstates of the destruction operator. 
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a = 2-^'\(3Q + iP / (3) (17) 

for the reference harmonic oscillator: 

H = PV2m + mujl Q"^ /2 = u,. a^a + h/2 , (18) 

with P = {mcUr-Y^'^. They can be obtained by displacing the normahzed ground state, 
|^x=o) = |0), of (|1^) to the phase space location X = {Q,P) according to, 

|f]x) = exp{{t/n){PQ - QP)}\0) . (19) 

Therefore, the Husimi function for pure states of 1-D system is, 

where we include the prefactor, for the case of normalizable states, to make the Husimi 
function integrable to unity over the whole phase space. 

The Husimi representation can also be viewed as a Gaussian smoothing of the Weyl 
representation [§, so the Husimi function is related to the Wigner function in the form, 

^^(X) = ^/dxW^(x)exp{-i||(x-X)||5} (21) 
where the "/3-metric" is defined as 

II (x - X)||, ^ [P\q - Qf + {p- Pf/n" ■ (22) 

This expression is the starting point for our approach developed in the next sections. Since 
WQy^ = l/vr^exp | — ^||(x — X)|||| is the Wigner function for a coherent state (see e.g. 0), 
(pip defines the positive definite projection, (0), of the Wigner function onto coherent states. 
Note that the Weyl representation is invariant under symplectic transformations (linear 
canonical transformations). The introduction of a metric in ( pTf ) implies that symplectic 
invariance does not carry over to the Husimi function. 

The fact that the coherent states are eigenstates of the destruction operator gives them 
analytical properties which are translated to the Husimi function We can separate 

the analytical part of the Husimi function if we use the unnormalized coherent states \z), 
such that l^x) = exp{—zz /2h}\z) , i.e., 



where the coordinates, 

z = 2-'/\[3Q-tP/(3) , z = 2-'/\[3Q + tP/(3) , (24) 

represent the complex phase space for 1-D systems. 

The function (zlip) corresponds to the wave function for the quantum state in a represen- 
tation of quantum mechanics introduced by Bargmann [|l^] (in the case of the Heisenberg- 
Weyl group), where the basis for the Hilbert space are the unnormalized coherent states 
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(althought overcomplete). These wave functions are holomorphic (entire) functions of the 
variable which behaves as a phase space coordinate. From the WKB construction in the 
Bargmann representation, Voros |^ derives a semiclassical approximation for the Husimi 
function in 1-D systems. This is presented in Appendix for the energy eigenstates of the 
problem of a particle under the action of a constant force and the results are compared with 
our approximation of Section 0. 



III. THE PARTICLE IN A BOX. 



To introduce our approach let us study the simplest case of a particle in the symmetric 
classical potential: 







i +00 





< 



V(g) 

with the pure states given by even eigenf unctions 
(2//)i/2cos(p,g//i) |g|< 

|g| > 



(25) 



(g|^n) 



with Pn = 'n'filn + 1) / 1 {n even) 



(26) 



Thus the semiclassical limit for a given classical momentum, pn, corresponds to the limit of 
large n. 

The Wigner function for the state is zero outside the box, whereas inside: 
Wjq. p) = ^\ ^'^f " + ^'^^ + P-^y^^^ + 2 cos i2p^q/n) ^^Py^^^ ] (27) 

21 [ T^[p-Pn) T^{P + Pn) T^P J 

where ?/ = | — (3'ifO<g<| and ?/ = | + g if — | < g < 0. 

The Husimi function can be calcutated in terms of the error function $(-2) (see Appendix 

0), 



2h 



Z2 



Z3 



2h 



Z4 



2h 



+ 2e 



2h 



Z2 



2Ti 



2h 



Z4 



2h 



(28) 



where the arguments of the error functions are measured from the "corners of the phase 
space box" shown in FIG.|, = /?(f + Q) - i(P - Pn)//3, z^ = /3(| - Q) + i{P - Pn)//3, 
^3 = /3(1 + Q) _i(p + p^)/^ and Z4 = (3{^-Q) +i{P + Pn)/f3; (the overlines in (|2S|) indicate 
complex conjugation). As opposed to the Wigner function, (p8|) is not zero outside the box. 
A semiclassical analysis of this expression can be made with the help of the asymptotic 
expansion of the error function $(iy) for \w\ sufficiently large, i.e. h ^ (see Appendix 0). 
If we replace each $(w), by the first term of this expansion. 
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f 1 - ^ ^e(w) > 1 
1-1-^ 5ReH«-l 

we get a good approximation of (pSf) except in narrow margins along the lines Q = —1/2 and 
Q = 1/2, that contract when h —>■ 0. This can be observed, for the region of interest inside 
the box and between the branches of the classical trajectory by comparing the plots (a) and 
(b) in FIG.|^. Furthermore, the function approaches unity, in the limit \w\ +oo, in 

the region |53m(iy)| < ^e{w) {^e{w) > 0). The intersection of these regions for each error 
function in (^) defines a central rectangle shown in the plots (a) and (b) of FIG.^. Hence, 
in the semiclassical limit the Husimi function is well represented by, 

1 f (P-Pn)^ (P+Pnf (f^+Pn) "I 

H{Q, P) ^ = <^ e + e + 2 e cos{2pnQ/h) } = 



2lpV7rh 



1 _ 
=e «^ 



1(3^71% 



e'lp |cosh(2p„P//32^) + cos(2p„Q/;i)} , (30) 



within that rectangle. This expression is explicitily positive everywhere except in its zeroes 
that lie on the axis P = where it simplifies to 

H{Q,P = 0) ^ -^e-^ {cos\pr,Q/h)} . (31) 

Alternatively, since the Husimi function is a Gaussian smearing of the Wigner func- 
tion (pID, it is possible to obtain a semiclassical approximation of (|28|) by performing 
the Gaussian smoothing over an expression that mimics the behavior of the Wigner func- 
tion in the semiclassical limit. With the help of the limiting form of the delta funtion, 
5{x — Xq) = lim^^^+oo sin[L(x — Xo)]/'/r(x — Xq), we observe that, semiclassically, the skeleton 
of the Wigner function (pT]) is 



W{q, p) ~ ^ {5 (p - Pn) + S{p + Pn) + 2 S{p) cos{2pnq/h) } . (32) 

Replacing this expression in (^) and limitting the integration to the range (— //2, //2), leads 
to 



Hscl{Q.P) 




(P+Pn)'^ \ 




nffi 










(33) 



where z[ = /3(| + Q), 4 = /5(f - Q), 4 = /3(| + Q) - ipn/P and 4 = /5(| - Q) + ipn/p. 
In this case we can replace each error function by unity only in the central region (between 
the vertical lines through Qi = —1/2 + p„//3^ and Qu = 1/2 —pn/ (3"^ in FIG.|(c)) and hence 
recover (0). Therefore, we also recover the position of the /c'th zero along the axis, 

Qfc = (2fc + 1)^, (34) 
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according to (|3TD, provided that Qi < Qk < Qu- In FIG.|^ we compare, on the Q-axis, the 
numerical computation of the Husimi function (|28|) , the asymptotic approximation based on 
(p8|) with ( p9D and our approximation based on the simphfied Wigner function (pS]). 

The form of (|30D indicates why the Husimi zeroes are hneary distributed inside the energy 



curve (in this case, the phase space box with \pn\ = V^mE). Indeed, the term with the 
hyperbohc cosine takes its lowest value along the Q-axis, which coincides with the amplitude 
of the oscillatory cosine term. Away from P = 0, the hyperbolic cosine dominates the sum, 
descending to a valley along this axis. The zeroes along the valley are determined by the 
minimum value of the cosine. The valley is very shallow because the Husimi function decays 
exponentially away from the classical region, but we can still evaluate its local minima. The 
order for the spacing of zeroes is 0{h), given by the phase of the cosine term. 

As autonomous 1-D systems always have integrable classical dynamics, the zeroes of the 
Husimi functions lie over lines, as suggested in That these lines, inside the energy 



curve, are valleys of the Husimi function is a general characteristic of these systems, as we 
will see in the following sections. It also seems to be a general characteristic that, when 
the number of zeroes is great, these valleys bifurcate for bounded states in systems where 
the curve of constant energy is closed. For fixed energy, in systems with bounded states, 
it is expected that the semiclassical approximations works well for large quantum numbers. 
As the number of zeroes of the Husimi function grows with the quantum number 0, these 
bifurcations should typically appear in the semiclassical regime. These bifurcations seem 
also to appear close to the energy curve. 

For the box, the fa//e?/ bifurcates close to the points Qi and Qu of the Q-axis (FIG.|]). Our 
approximation ( |33D describes these bifurcating valleys , althought without any oscillations 
to indicate the presence of zeroes (see plot (c) of FIG.^. The absence of zeroes in these 
valleys shows that the expression ( P^D for the Wigner function close to the edges of the box 
is not valid. In fact, the Wigner function (^7]) decreases to zero close to the edges of the 
box and is strictly zero over them. In constrast, the expression ( P^ does not decrease in 
the direction of the Q-axis and does not vanish over the edges. 



We stress that the only approximation used to obtain (|33|) is to take ( p2|) as the Wigner 
function. Moreover, by extending the limits of integration to infinity, in the smoothing 
of the Wigner function, we obtain the expression ( ^Of ) for all points inside the box. This 
approximation has no bifurcating valleys at all. But, as we knew, this approximation is not 
valid close to the edges, because making the limits of integration go to infinity is equivalent 
to making the error functions in (^) approach unity, which is not valid in the region where 
the bifurcations occur. 



IV. GEOMETRICAL APPROACH. 

In the last example we obtained an approximation to the Husimi function by perform- 
ing the Gaussian smoothing over an expression that represents the skeleton of the Wigner 
function in the semiclassical limit (0). It was shown that this provides the general behav- 
ior inside the energy curve and allows us to obtain the distribution of the Husimi zeroes, 
although not close to the edges of the box. Here we implement a similar approach for the 
Husimi function of energy eigenstates in systems where Berry's semiclassical approximations 
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for the Wigner function are valid (see Section H). 

The ideal semiclassical approximation to the Wigner function used in the smoothing 
should be Berry's uniform approximation, that represents the oscillations inside the energy 
curve S, and is uniformily valid along it. However, the integration would be very difficult 
to handle. The skeleton of the uniform approximation essentially consists of an Airy peak 
close to the curve S, that in the classical limit turns out to be a delta function ([Tl|) along 
it, and oscillations inside that are well represented for the Berry's simple approximation 
away from S. Semiclassically, as far as integration is concerned the Airy peak is equivalent 
to the delta function. So, for evaluation points, X, close to the energy curve, the Husimi 
function is well represented by the integral, 

/(X) = ^-^^i— ^ /dx<5[H(x)-E] exp{-i||(x-X)||j} . (35) 



This integral is everywhere positive, smooth, and peaked along the energy curve S. Evi- 
dently, this integral is dominated by the region where ||(x — X)||^ is a minimum: approxi- 
mately ^^7= exp | — i II (xc(/3) — X)||^|, where Xc(/3) is the point on S closest to the point X 
in the sense of the norm ||(. . .)||^. Thus, close to the energy curve, this integral is essentially 
the Gaussian semiclassical approximation around the torus, first encountered by Takahashi 
1^ in a geometrical approach, and rederived by Kurchan et. al. [ffH] in the context of 



the Bargmann representation. These approximations have no oscillations to indicate the 
presence of zeroes. So, placing the evaluation point far from S, we add to the integral of 
the delta function, a local integral over the simple approximation. This is in the spirit of 
expression (^) for the problem of a particle in a box. Indeed, the only difference is that 
the cosine oscillations are now spread within the interior instead of concentrated as a delta 
function along the g-axis because of the particular torus geometry. 



The integration (|2T| ) over the simple approximation, inside £, can be performed ana- 
lytically by making some further approximations (for the details see Appendix P). The 
Gaussian function in (Ell) defines an effective area for the integration centered on X. This 



effective area can be characterized as the area of value 2nh, enclosed by the ellipse, 

i||(x-X)||J = 2 (36) 

So, inside this area we can approximate the action in (|D by 

^(x) ^ 5(X) + e(X) A (x - X) , (37) 

in the semiclassical limit. Since the denominator in (^) does not depend on h, we take the 
simplest approximation, 

D(x) ^ D(X) . (38) 
Then, the result for our approximation to the Husimi function is, 

Hscl(-X) = 7= , COS < — 7 > + i X . 39 
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As we found in the last section, the Husimi zeroes inside the energy curve are located on 
a valley for 1-D systems. The approximation (p9|) contains all the geometrical ingredients 
to understand the origin of this valley. This expression generally has minima rather than 
zeroes in the valley. These minima approach the Husimi zeroes when h ^ 0. However, 
some of these minima could be negative in this approximation. This problem can be fixed 
if we include the second order approximation, ^ (x — X) 7i (X)(x — X)*, in the expansion 
for center action in (|57D, where the Hessian matrix is 



9pqS - 



dq{Q 



dgpS 

dppS = 



(40) 



(we have applied the relations (|T5D for the gradient of S{x.)) and t denotes the transpose. 
Hence, our refined approximation 



HscL (X) 



exp{-e(X)/4^} 



UJ 

\siiq 



X 



X COS 



D(X)|det^(X)| 

arg[det^(X)] j ^ ^^^^ 



TT $(X 

4 ^ 



(41) 



Here, .4.(X) is the complex matrix, 

^(X) = 





1/13^ 



+ 2^(X) 



(42) 



the argument of the exponential is 

_ lle(X) 11^(1- det n{X)lA) + \ e(X) n{X) e*(X) {d,{Ql2(3^ - P'd,{Ql2) 
- |det^(X)P 

and the phase in the cosine is 

-||e(X)||^ {d,{Q/2(3'-(3%{Q/2) + la^)nmem (l-det7^(X)/4) 



(43) 



$ X 



det^(X)|^ 



(44) 



If we only use the approximation to first order for the center action in (^), the Hessian 
matrix is zero, so, det^(X) = —1, 6(X) = ||^(X)||^ and $(X) = and we recover the 
approximation (0). For the cases of non-normalizable states the prefactors in (pSj), ( pOf ) 
and (|4l|), change according to the definition of the normalization of this type of states (see 
Section |T| and Section |V| for an example). 

The expressions (^) and (^) are valid only inside the energy curve and depend only 
on the properties of the curve S, like the semiclassical Wigner function. The second order 
approximation to the center action, that yields our approximation (^T]), only provides small 
corrections to the argument of the exponential and specially to the phase in the cosine. The 
corrections of the phase in the cosine improve the position of the minima over the valley and 
hence the approximation to the zeroes. 



12 



The fact that the denominator in ( p9D vanishes on the curve £ is not a major problem, 
because only the tips of the valley are close to the energy curve, where the simple approx- 
imation plus the delta function ( ]Tl| ) is not a good representation of the behavior of the 
Wigner function. So, our approximations (|39|) and (^11) do not hold close to the curve £. 
The evaluation points x, of the Wigner function, that effectivelly contribute to the smooth- 
ing (p^, are enclosed by the ellipse (|36D centered on the evaluation point, X, of the Husimi 
function. We predict that the approximation based on the mimic Wigner function breaks 
down wherever the ellipse enclosing X approaches £. The shape of the ellipse depends on 
the Husimi parameter /3, so that the region of validity of the geometrical approximation will 
be parameter dependent. 

We now discuss the fact that the geometrical approximations (^) and (|4l|) contain the 
contribution of a single chord, even though the points within the Wigner caustic C are the 
centers of three chords. This is simply due to the Gaussian dependence on the chord length, 
defined in (^^, which allows us to keep only the shortest chord in the Husimi function. 
Furthermore, we need not consider the caustic itself, because, as we will see, the valley 
of zeroes is not affected by it. Even though the simple approximation breaks down along 
it, by predicting a spurious singularity, the correct finite Airy peak along this line does 
not counterbalance the fact that the coalescing chords responsible for this catastrophe are 
longer than the normal third chord. This is because the curve £ is a locus of maximal chords. 
Therefore the Husimi function will be dominated by the single normal chord on £, which 
generates cosine oscilations well described by the simple theory. Hence, the sum over the 
differents chords, centered on points on the Wigner 's caustic C and inside it, that appears 
in the simple approximation (^, is not necessary in (^) and (|^). We only need to consider 
the chord that is continuous through each of the two sides of C that are crossing the valley. 
For this chord the denominator -D(X) does not diverge (see FIG.[I|). 

In the next sections we present two examples to show how (^) and (^9l) operate. The 
first example is an unbounded problem whose convex energy curve is open and without a 
Wigner caustic. In this problem most of the calculations can be made analytically. The 
second is a bounded problem with a closed energy curve that is smooth and convex. This 
example has a Wigner caustic. In this case all the calculations were numerically. 



V. PARTICLE SUBJECT TO A CONSTANT FORCE. 

Let us apply the approach described in the last section to the problem of a particle 
under the action of a constant force, F, that is to say with the classical Hamiltonian H(x) = 
p^ /2m — Fq. This is an unbounded problem with continuous energy spectrum where the 
eigenfunctions can be normalized to a delta function in E {i.e., J {'ipE'\^){^\'^E)dq = S{E — 

E)) i, 



{(i\i^E) 



|F|V2 



2m\F\ 



1/3 



Ai<-(g-gr 



2mF 



1/3 ~ 



(45) 



Qr = —E/F is the turning point of the classical trajectory for an energy E. The Wigner 
function is in this case Il5|,|6|l 
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Vr(x) 



5m 



1/3 



Ai 



em 



nl/3 



(H(x) - E) 



(46) 



It is easy to see that this expression coincides with (H) (the prefactor equal to unity in this 
unbounded problem, according to the normalization choosen above). So, the Transitional 
approximation to the Wigner function is exact in this case. 

The Husimi function for this problem can be caculated analytically (Appendix |C|), the 
result is. 




(47) 



where = qr + TnF/2f3'^ and |-Bp is the normalization constant. The distribution of zeroes, 
in the concave side of the curve S, of constant energy, is shown in FIG.^. Since the zeroes 
are those of the Airy function in (^Tj), which only occur for a negative real argument, their 
distribution is along the Q-axis for any h value. For an energy E the classical turning 
point is fixed and so is Q^. Thus, we only have to change the scale for the X coordinates 
to make the Airy function in (^) invariant with h. Due to this scaling property, we can 
analyze semiclassically the distribution of zeroes by fixing h and by looking at the behavior 
of the Airy function for values of Q far away from Q^, on the concave side of In this 
region we can replace the Airy function in (^) by its asymptotic form (|^), where now the 
argument is complex. The result is the same if we construct the Husimi function in (p3D, 
with the approximation to the Bargmann function (plO|) , obtained in Appendix by the 
complex WKB method. Therefore, the direct semiclassical analysis of expression {^7\), or the 
application of the complex WKB method, bring about the same distribution of the zeroes 
in the concave side of S. This distribution is given by the zeroes of the cosine in ( |D10| ) over 
the real axis, 



[37rM2fc + 3/2)]^/3 



k = 0,1,2, 



(48) 



where now F > 0. The order for the spacing between the fc'th and the {k + l)'th zero is 
0{h'^^^). However, these zeroes accumulate on as h 0, so it is more relevant to derive 
the asymptotic spacing near a fixed position Q. Approximating the Airy function in (^71) 
by its large argument from (|^), we then obtain the spacing of minima as 0{h) in agreement 
with Leboeuf and Voros p!6| . 

To compare these results with the geometrical approximation, we note that for the par- 
ticle subject to a constant force, the centre action is 



S{X) 



1/2 



-iH{X)-E)] 



3/2 



(49) 



and the denominator in 
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D(X) = - = P— [-(HiX) - E)Y^' . (50) 
m \ m 



jl/2 

Hence, in our approximations (PU[) and (HTf) the first term can be calculated analytically, 



where now the prefactor is (S/vr/i)^/^. The classical limit (0) can be obtained by applying 
formula (0) to the Wigner function (^61). The integral (^) over £ becomes, 

/(X) = ^ / rfpexp [-\mqE{p) - Qf + {p- P)l(3']] (51) 

where qE{p) = p^/2mF + g^- This is a non-oscillatory smooth function, peaked on S, that 
decreases monotonically away from the energy curve. 

The geometrical origin of the valley of zeroes along the axis P = 0, in the concave side 
of S, can now be understood with the help of our approximation (pOl) . In fact, from the 
argument of the exponential, we see that for values of (3 allowed by our approximations (see 
Section |rv| ), ||^(X)||^ is essentially equal to the square of the chord's length. Along the 
Q-axis, FIG.^ shows that the prefactor of the cosine in (|39|) has almost the same value as 
the integral (|51|). Away from this axis, the length of the chord grows, making the oscillatory 
term so small that the second term dominates the sum, creating in this way a valley. Along 
this valley the oscillations of the cosine generate a sequence of local minima. However, the 
position of these local minima are shifted, relative to the position of the Husimi zeroes, 
approximately by the distance (Q* — Qr) (FIG.^. Moreover, for points Q far away from 
Q*, these local minima become negative because the prefactor of the cosine becomes greater 
than the integral (see FIG.^. 

The corrections given by our second order approximation ( pj ) fix these problems. In 
fact, the corrections to the argument of the exponential, given by (^Hf) , ensure that the 
local minima are positive on the axis and the corrections to the phase in the cosine, given 
by (|4^) , improve the position of the minima relative to the Husimi zeroes (see FIG.^ and 
FIG.|^). We also compare, in FIG.^, the general behavior of ( |^ ) and the Husimi function 
calculated numerically on P = 0. One observes a general agreement that improves as Q 
recedes from Q^, (equivalent to the limit ^ ^ 0, in this problem). In this semiclassical 
limit, the same figure shows that the minima become zeroes. FIG.^ shows the relative error 
between the position of the Husimi zeroes (calculated numerically), the position being given 
by the minima of ( PD and (|39|) (shifted by the distance {Q^ — q^)) and the zeroes (^8]). As 



expected, our approximation (^) does not work well close to the energy curve where the 
mimic of the Wigner function used for the smoothing is not a good approximation. 



VI. GENERIC CASE. 

In this section we apply our geometrical approach to the problem of a particle subject 
to an asymmetric anharmonic potential. The classical Hamiltonian is 

This system is an example of a general system with a convex, closed energy curve £ and a 
Wigner caustic C FIG.|1| shows typical curves, 8 and in this system. 



15 



In this bounded problem, we fixed the classical energy curve at i? ~ 30.8175 and we 
calculated numerically the distribution of the zeroes of the Husimi functions inside it for 
two energy eigenstates. The latter correspond to two values of quantized h: the eigenstate 
n = 30 for a value of ^ ~ 0.508236, and the second is the eigenstate n = 45 for a value of 
h ^ 0.340691. In order to simplify the calculations, all the others parameters of the problem 
(including /5) were set to unity, except go = 4.0 and A = 0.1. 

FIG.^ shows the distribution of the Husimi zeroes for these states. The zeroes are dis- 
tributed along lines, as expected for a system with integrable classical dynamics |]16| . These 



lines are very shallow valleys of the Husimi function. Since the energy curve is symmetric 
with respect to the Q-axis, the distribution of zeroes also mantains this characteristic. In 
Section |T| we anticipated, as a general charcteristic, the bifurcation of the valleys in the 
semiclassical regime for bounded states in systems with a closed energy curve. Here, we 
have a generic example where a principal valley bifurcates in each half plane (FIG.^. The 
asymmetry in the lengths of these bifurcating valleys reflects the asymmetry of the curve S 
with respect to the P-axis. For each quantum number n, the majority of the Husimi zeroes 
belong to the principal valley. For the parameter j3 chosen, the principal valley runs parallel 
to the vertical side of the caustic. At the middle it passes very close on the outside of C and 
them crosses the cusps of the caustic (see FIG.^. This shows that the distribution of the 
zeroes is not affected by the presence of the Wigner caustic. 

Our approximation (pOf ) to the Husimi function inside S, supplies the geometrical insight 
for the origin of the principal valley. As we saw in the example of the last section, this valley 
is located where both terms of the approximation are of the same order. The amplitude of 
the oscillatory term is dominated by the exponential, whose argument is essentially equal 
to the square of the chords length. If we follow the chords length along the level curves of 
the center action, that is, essentially the phase curves of the cosine in the oscillatory term, 
we observe local minima of this length restricted to these curves. Away from the valley the 
length of the chords grows, making the oscillatory term so small that the smooth second term 
dominates the sum. This is also what happend in the problem of the constant force (Section 
^ where the Q-axis is the locus of minima of the chords length restricted to the level curves 
of the center action which cross the axis orthogonally. Along the valley, the oscillations of 
the cosine produce a serie of local minima of (^) that indicate, in a first approximation, 
the position of the Husimi zeroes. The position of these minima are very close to the points 
where the cosine takes its minimum value, only slightly modified when we consider the sum 
of the two terms. In FIG.^ (a) and FIG.|TD| (a) we illustrate, for each quantum number, 
the geometrical method to locate the valley and the Husimi minima. We display the level 
curves of the center action for which the cosine takes its lowest value, and the level curves of 
the chords length. The tangency of the two sets of curves determine the restricted minima 
of the chords length along the chosen level curves of the center action. These points are 
candidates to be, approximately, the local minima of ( p9D after performing the sum of the 
two terms. The valley of minima passes through all these points in this approximation. 
The comparison in FIG.^ between the valley of minima of (PD|) and the zeroes of the Husimi 
function shows that our approximation to the principal valley of zeroes holds farely well until 
the bifurcation. Although this approximation to the valley continues after the bifurcation, 
it does not take into account the bifurcation itself. Even this continuation, is no longer such 
a good approximation of the longer valley beyond the bifurcation. 
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The local minima of (^) along the principal valley have almost the same spacing as the 
Husimi zeroes. However, the absolute position of the predicted zeroes is not so precise (see 
FIG.|^ (c) and FIG.|10| (c)). Furthermore, some of these minima become negative, because 
the prefactor in the cosine becomes greater than the integral in the second term of the 
approximation. We found the same situation when we applied (pOf) in the last section, so we 



again make use of the refined approximation (0). This supplies corrections to the chord's 
length in the argument of the exponential and to the center action in the phase of the cosine. 
Therefore, we can use the same geometrical method to find the approximate position of the 
local minima of (51) of the Husimi zeroes in the semiclassical regime. Hence, FIG.^ (b) and 
FIG.|TO| (b) display the phase curves of the cosine for minimum values, and the level curves 
of the argument of the exponential. The tangencies of these two sets of curves determine 
the restricted minima of the argument of the exponential in (p]), along the level curves of 
the phase of the cosine. Clearly, these points belong to a valley because, away from the line 
that passes through all the restricted minima, the smooth second term in (5T) dominates the 
sum exponentially. Moreover, since the cosine takes its lowest value at these points, they are 
close to the local minima of (|^). The comparison of these points with the Husimi zeroes in 
FIG.|^ (d) and FIG.|10| (d), shows that they are a good approximation to the zeroes along 
the principal valley until the bifurcation. 

Our approximation (^T]) also fails to take account of the bifurcation. The valley of local 
minima is again a very good approximation to the principal valley of the Husimi function, 
but now, it also represents accurately its continuation along the longest bifurcating valley 
(see FIG.^. However, along this bifurcating valley there are no local minima of (0) to 
indicate the presence of zeroes, because there the prefactor of the oscillatory term is much 
smaller than the integral of the second term. This can be observed in FIG. 11 and FIG. 12 
where we display the logarithm of our approximation and the Husimi function along the 
valley of local minima of (^if ), for each quantum number. 

We saw that the valley of zeroes of the Husimi function is not affected by the Wigner 
caustic. The principal valley, that runs parallel to the vertical side of the caustic on the 
outside of C , crosses the cusps of the caustic at the tips of this side. This is reflected in our 
approximations (^) and ( ^T]) since only the contribution of a single chord inside the Wigner 
caustic is needed. Because of the dependence on the chord's length in our approximations, 
we showed that we only need to consider the shortest chord that is continuous through £, 
irrespective of the two chords that coalesce along this caustic. This shortest chord generates 
the valley that runs parallel and very close to the principal valley of the Husimi function 
unaffected by the close lying Wigner caustic (see FIG.H). 

To end this section, we note that the consideration of states corresponding to a flxed 
energy, that are quantized by varying h, leads to a spacing of the Husimi zeroes of 0{h). 
This is easily seen in our simple approximation (|39D , for a flxed location along the valley 
that is classically determined. The wavelength of the oscillations that determine the minima 
is proportional to h, whereas ^(x) is a classical action. 
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VII. CONCLUSIONS. 



The semiclassical approximation of the Husimi function has been derived by integrating 
the semiclassical Wigner function with a Gaussian window. This is not fundamentally 
different from the calculation of probability densities of position or momenta as projections 
of the Wigner function, except that now we project onto coherent states. In each projection, 
we obtain a classical approximation by substituting the Wigner function by a delta-function 
along the classically allowed region. Here, this leads to a narrow ridge along the classical 
region, which is supplemented by an oscillatory term derived from the centre and chord 
structure within the energy curve. The oscillations of the latter along the classical shallow 
valleys combine to form local minima, that indicate the positions of the Husimi zeroes 
inside the energy curve. This geometrical explanation of the linear distribution of minima 
cannot be pushed to the point of predicting absolute zeroes, but it may be nonetheless 
surprising that their positions are asymptotically accurate, though obtained by subtracting 
two exponentially small terms. 

The advantage of deriving the intermediate approximation (|39|) is that the location of 
the valley has a simple dependence on the minimal chord along curves of constant centre 
action S'(x). This curve is purely classical, once the excentricity j3 of the coherent states 
defines the phase space metric. The corrections added to our complete formula (^Tj) are also 
classical. They mostly alter the distribution of zeroes along the valley, so we find that the 
valleys are basically determined by the classical structure, in agreement with Leboeuf and 
Voros [|16] . 



It is perhaps surprising that the Wigner caustic C does not affect the position of the 
Husimi zeroes. However, this fact is in agreement with previous calculations for the projec- 
tions of the Wigner function . In each case the integration singles out a contributing chord 
from the semiclassical Wigner function, while ignoring the other possibly singular chords. 
Thus we can understand the complexity of the Wigner function as arising from the necessity 
to account for diverse square integrable projections. 

We have limitted our considerations to autonomous Hamiltonian systems with one de- 
gree of freedom, which are necessarily integrable. Our approximations are equally valid for 
integrable classical maps on the plane onto itself: the zeroes are always predicted to lie along 
the locus of minimal chords, the minimum being evaluated along lines of constant phase for 
the Wigner function. 

The chord structure generalizes to tori of higher dimension 0] , so our methods will also 
be extendable to the study of Husimi functions of integrable systems with more than a single 
degree of freedom. In particular they may help to define the zero-manifolds in this case. 
For chaotic systems, we know that the chord structure is also present, though it involves 
individual orbits [§]. The challenge lies open, to explore the relation between the structure 
of the Husimi and the Wigner functions for nonintegrable systems. 
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APPENDIX A: 



The Husimi function for the problem of a particle in a box with hard walls. 



Here we show the fundamental steps in the derivation of (pS]). We start with the expres- 
sion of the Husimi function given by the formula (^OD , 



2nh 



+00 



(^x|g)(g|^n)rfg 



li-Khf/^ 



, (Al) 



where {q\Qx.) is the normalized coherent states in the position representation (see for example 
P^ ) and (qli'n) the even eigenfunction (PB]). If we express the cosine in the last integral as 
{l/2){e'P^i/^ + e-*P"^/^) we have 



P 



4i{7ih)y^ 



e 2ftf^ e 



.■(f-Pn)Q 



+1/2 

I 

4/2 



-\^^W(q-Q)+i{P~Pn)/l3]y 



dq + 



_(P+Pn) . (P+p- 

+ e 2fi/32 e-« 



i-i/2 



(A2) 



Making a change of variables in the expression between braces in each integral and using 
the definition of the error function 



TT JO 



'dy 



the Husimi function becomes 



HiX) 



1 



UpVTrh 



2/1 



2/1 



2/1 



(A3) 



(A4) 



With the help of the identity |wi +tf2| 



|w2p+ 2^e{wiW2) of the complex numbers 



we arrive to the expression ( pSD for the Husimi function in this problem. 

Finally, we give the asymptotic expansion of the error function $(w), used in Section 
PI, in view of the confused and even incomplete form that it appears in the usual references 
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$(w) 



-1 - 



[Fn{w) + 0(|«;|-2("+i))] ^e{w) > 
[Fn{w) + 0(|w|-2("+i))] sRe(u;) < 



\w\ large 



(A5) 



where Fn{w) = ELo 



(-l)fe(2fc-l)!! 
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APPENDIX B: 



Details of the geometrical approximation to the Husimi function within the energy curve. 

Here we derive the oscillatory term of our expression (0). We start with the Gaussian 
smoothing (^) over the simple approximation (H) within the energy curve. If we apply the 
approximation ( pT]) to the center action, induing the second order term ^ (x — X) 7i (X)(x — 
X)* (with Ti (X) the Hessian matrix (^OD), and the approximation (|38D for the denominator, 
we have 



7rv/2^(cu-i)7r;i^D(X) 
where X(X) is the integral. 



exp 



h 



4 



J(X) 



(Bl) 



J(X) = / dq f dp eM-"i('?-Q)2-a2(p-P)2+a3{g-Q)(p-P)+a4(g-Q)+a5(p-P)] ^ (^g2) 
J — oo J—oo 

with the complex coefficients: ai(X) = (3"^ — idg{C,p)/2, a2(X) = 1//?^ + idp{^q)/2, a3(X) = 
i dp{^p), a4(X) = i and a^iX.) = —i ^q. This double Gaussian integral can be performed: 



X(X) 



ai{a2 — a|/4ai)| 



,2 

'A 
ai 



(02 - ai/4ai) 



+ 02) 



(B3) 



where 6^1 = arg(ai) and 62 = arg(a2 — 03/401) with tt/2 < 9i, 62 < 7r/2. This result can be 
written in a more elegant way with the help of the complex matrix, ^(X), defined in (ff3). 



J(X) 



det^(X)| 



exp 



1 e(X) ^(X) e*(X) . arg(det ^(X)) 1 



Ah det ^(X) 



1 



(B4) 



Replacing (P4|) in (|B1|) and solving the real part we find the oscillatoty term in (P|). 



APPENDIX C: 

The Husimi function for the problem of a particle subject to a constant force. 

The idea is to solve the Schrodinger equation for the eigenstate \iPe)-, written in the 
Bargmann representation, to find the Bargmann function in (|23| ) for this problem. To write 



the equation in the Bargmann representation, we proceed in the standard form: we write 
the quantum Hamiltonian as a function of the creation and destruction operators ([T^) in 
normal order [i.e. the operators to the right of the d operators) and then using, 

{z\d\i,) = hdMi^) {z\d^\ij) = z{z\i;) , (CI) 

(where the are the unnormahzed coherent states defined in Section |I|) we get 
f /?2 F Fh ^ 
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This equation can be put in the form 

{[d,-f{z)]^ + ciz-z,)}{z\ijE) = 



(C3) 



where c = {l/h^)8mF/ /3^^/2, z^ = {/3/V2){qr.+mF/2(3^) (with = -E/F the turning point 
of the classical trajectory of energy E), and the function f{z) = {l/h){z- {2mF / (3^ ^/2)) . 
The general solution of ( PSf ) is 

(C4) 



= exp{(7(z)} Ai{-(^ - z*)[c]i/3| ^.^^^ Q^g^^-^ ^ f^^-^ 



Therefore, the Bargmann function for an eigenstate of the problem of a particle subject to 
a constant force is 



{z\iIje) = Bexp 



1 z' 



2mF 



n\2 ^/2f3'-- 



z\ } Ai < — (z — z* 



8mF 
¥V2i3\ 



1/3^ 



(C5) 



where i? is a complex constant. Replacing (|C5|) in (E3|) we obtain the Husimi function (B^. 



APPENDIX D: 

A semiclassical approximation to the Husimi function for the problem of a particle subject 
to a constant force through the WKB method in the Bargmann representation. 

Here we follow the WKB construction in the Bargmann representation given by Voros ^ 
to obtain a semiclassical approximation to the Bargmann function {z\iPe) and then, through 
(^), a semiclassical approximation to the Husimi function. 

We use the WKB construction based on the Weyl symbol i/i^(x), of the quantum Hamil- 
tonian H, for the cases where it does not depend on ^ {i.e that it coincides with the classical 
Hamiltonian, H). It is easy to verify, through (^, that Hw = H when the classical Hamil- 
tonian is of the form @. As this is the case for this problem, we write H instead of Hw- 

Voros argues that to build a WKB solution of the equation B.\iPe) = EIi/je) in the 
Bargmann representation, we should just apply the same formulae as for the Schrodinger 
representation [i.e. the common position representation), while replacing x = {q,p) — > [z, z) 
and h ih. Therefore, the eigenvalue equation admits local asymptotic solutions to leading 
order in h, 

-1/2 

exp{S{z)/h} , (Dl) 

where ze{z) is the function defined implicitily by the relation that defines the classical energy 
curve in the [z, z) coordinates, 

H(z,z)=E, (D2) 
and S{z) is the classical action in the complex coordinates. 



We) 



9H 



dz 
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S{z) 



ze{z )dz . 



(D3) 



The fact that z and its complex conjugate z are related themself by H(z, z) = E restricts z 
to the real energy curve. This means that the function ze{z) is a branch of ( p2| ), defined in 
a sheet to which the energy curve belongs, and is single valued over it. For nonanalytic H, 
there is no guarantee to have an analytical continuation of ze{z) anywhere outside. Thus 



(Dl) is well defined only for z on the real energy curve and it is globally regular, since there 
is no turning point anywhere. 

In order to obtain a holomorphic approximation far away from the energy curve we use 
the fact that, in this problem, the classical Hamiltonian is analytic in both variables z and 
z, so the energy relation (^) defines implicitily ze as a multiply valued function of z {i.e its 
domain are the sheets of some Riemannian surface). However, outside the real energy curve 
ze{z) is no longer the complex conjugate of z, so we use a less confusing notation denoting 
as y the independent complex variable canonically conjugate to z, 



/3q 



■ P 



/3q 



■ P 



where q and p are now complex. Then, the complex energy curve in this problem is 

H(z, y) = aoy^ + ai{z)y + a2(^) = , 



(D4) 



(D5) 



with the coefficients: oq = —(3"^ /Am, ai{z) = f3'^z/2m — F/-\/2/5 and a2{z) = —(3'^z'^/Am — 
Fz 1^2(5 — E. This is an equation of degree two, so the explicit branches y = yE{z) are 
defined over a two-sheet Riemannian surface. If we make the simple change of variable 
w = 2aoy + ai, we obtain the equivalent equation 



w —uiz) 







(D6) 



where u{z) = of — 4aoa2 = —[V2i3F/m]{z — z^) {z^ = [f3/y/2]{mF/2(3'^ + g^) and = 
—E/F is the turning point of the classical trajectory of energy E). The function w{z) is 
defined over a Riemann surface with branch points z = z^, and z = +oo. The branches are 
Wj{z) = ■i[V2/3F/mY/^y/z- z^ and w„(z) = -i[V2l3F/mY/^^z - z,, {F > 0). Here we use 
the notation: ^/z — z^ = ^/re^^^'^, < 9 < 2tt; so we can also consider z in the ordinary 
complex plane and Wi{z) and Wii{z) as two different functions. Since y is a single valued 
function of w, the branches yE^z), or equivalently the solutions of (P5|) are, 

1/2 

^/Z — Z^ + Z — ^^i^T 

(F > 0) . 



yu{z) 
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1/2 



y/z - Z^ + Z 



2mF 
2mF 



(D7) 



The WKB approximation is then a linear combination of solutions of the type 
each branch, valid away from the energy curve: 



OR 

dy 
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dy 
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e-xY>{Sn{z)/h] 
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(D8) 
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with the complex actions for each branch 



yi{z)dz 



yn{z)dz 



(D9) 



Hence, the semiclassical approximation to the Bargmann function in this problem is 

{z\%1)e) eicpl-y — 
{ 2 

X cos 





'8mF' 








(DIO) 



This expression can also be obtained applying the asymptotic form (|^ to the Airy function in 
the Bargmann function (|C5| ) of Appendix |C[ Replacing ( plOj ) in (|23|) leads to a semiclassical 
approximation to the Husimi function. 

Since the zeroes of the Husimi function are the same of those of the Bargmann function, 
the semiclassical distribution of the Husimi zeroes can be obtained from ( |D1U| ). However, 
besides the zeroes (^) over the real axis, (piO| ) predicts spurious zeroes over the straight 
lines that start in z = and have the directions 6 = 27r/3 and 6 = Ait/ 3. Notwithstanding, 
considering the region of validity for applying (|^ in (|C5|) , we see that those zeroes are in a 
region where (|D1CI|) is not an approximation of the Bargmann function (|C5|) . 
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CAPTIONS 

FIG.l. Semiclassical geometry of the Wigner function in 1-D systems for a typical 
smooth convex curve, of constant energy. The full chord near the Wigner caustic £, 
corresponds to a center point outside it. As the center moves through £, a bifurcation 
occurs. When the center lies on C (dashed chords), a second chord is born. Notice that the 
phase space velocities at the tips of this new longest chord are parallel, cancelling the area 
(^), while for the shortest chord this does not happend. Finally, when the center is inside 
£, there are three chords (dotted). The other elements of the geometry are explained in 
the text. 

FIG. 2. Husimi function plots (on a logarithmic density scale) for an even eigenstate in 
the problem of a particle in a box with hard walls; the stress is in dark for the greatest value 
and in white for the lowest one. (a) Calculated numerically from (|20| ) with the function 
(f2x|V') expressed in the position representation, (b) Asymptotic approximation for the 
error functions in ( ^5]) given by (PPj). (c) Smoothing of simplified Wigner function, (P^. 
The center of the white spots represent the zeroes. The horizontal white lines are the 
branches of the classical trajectory = 0.3 between the limits of the box. The points at 
the tips of this branches are the "corners of the phase space box" cited in the text. The 
central rectangle encloses the intersection of the region |Q'm(ti')| < 3fte(w) with 3fte(tu) > 0, 
for all the error functions in (pS]). The vertical lines in plot (c) enclose the intersections of 
the same region for the error functions in (^31). 

FIG. 3. The Husimi functions over the axis P = for the plots in FIG.^. The dotted 
curve corresponds to the plot (a), the dashed curve corresponds to the plot (b) and the 
full line is for the plot (c). The dashed vertical lines mark the limit of validity of the 
approximation (pl|). 

FIG. 4. The Husimi function, on a logarithmic density scale, of an eigenstate for the 
particle under the action of a constant force. The parabola is the curve of constant energy. 
The center of the white spots on the axis P = represent zeroes. The circle in the bottom 
left corner represent the curve (^) where the Gaussian smoothing (|21|) is significative. The 
symbol (x) in P = indicates the point Q^,. The values of the parameters used are: E = 10, 
F = 1, n = 2, /3 = 1 and m = 1. 

FIG. 5. Numerical camparison, along the Q-axis, of the prefactors of the cosine in (^If ) 
(dotted line) and in ( pQ] ) (dashed line) with the integral (|5T| ) (full line). Notice that when 
Q increases, the prefactor in (^) is grater than the integral (|5TD . The vertical dashed line 
is at position and the dotted one is at Q^k's. 

FIG. 6. The logarithm of our approximation ( PD (black curve) and of the Husimi 
function for the problem of the particle under the action of a constant force (gray curve), 
along the Q-axis inside the energy curve. The sharp inverted peaks indicate the position 
of the Husimi zeroes in the black curve, and the minima of (|4ll) in the gray curve. The 
vertical full lines stress the position of the Husimi zeroes. The vertical dashed line is at g^'s 
position and the dotted one is at (5^,'s. (o) indicates the Q-positions for the minima of our 
approximation (^) and (x) for the minima of (0). The second line of (o) are the minima 
of ( p9D shifted by the distance (Q* — qr). (o) are for the zeroes ( ^8]) of the semiclassical 
Husimi function obtained by the WKB method in Bargmann representation. 

FIG. 7. Percentage relative error, (5fc/Afc_i^fc)100% , between the positions of the Husimi 
zeroes (calculated numerically) and the positions given by the approximations to the zeroes. 
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Sk is the distance, on the Q-axis, between the fc-zero of the Husimi function and the position 
given by one of the approximations. A^.i ^ is the distance between the k — 1 and the k-zero 
of the Husimi function. The zeroes are counting from left to right, (o) are for the case of the 
zeroes ( ^8|) of the semiclassical Husimi function obtained by the WKB method in Bargmann 
representation, (x) for the minima of our approximation (^) and the (o) for those of our 
approximation (^) shifted by the distance (Q* — qr). 

FIG. 8. Distribution of the zeroes of the Husimi functions of two energy eigenstates inside 
the energy curve with E ^ 30.8175, for the problem of a particle subjet to an asymmetric 
anharmonic potencial (^2]). The symbol (x) indicates the position of the Husimi zeros, (a), 
for the eigenstate n = 30 for a value of h ^ 0.508236 and, (b), for the eigenstate n = 45 for 
a value of ^ ~ 0.340691. S is the energy curve and C the Wigner caustic. The dashed line 
represents the valley of the Husimi zeroes. The dotted line represents the valley of local 



minima of our approximation (|39[) and the full line, the minima of (|41|). The circle at the 
upper right corner of each figure represents the curve (PB]) for the range of the Gaussian 



smoothing (pl|). 

FIG. 9. Geometrical method for locating the position of the local minima, along the 
valley of approximations (|39| ) (a) and (^T]) (b). The value of h corresponds to the quantum 
number n = 30. (a): dotted lines indicate the level curves of the center action, while the 
level curves of the chords length are given by the full lines, (b): dotted lines indicate the 
level curves of the phase of the cosine in (|4l|) , while the level curves of the argument of the 
exponential are given by the full lines, (c) and (d): the positions of the Husimi zeroes for 
the eigenstate n = 30 (x), and the approximate position of the local minima ((o) for (p9|), 
(o) for 0). 

FIG. 10. Idem FIG.^ for a value of h corresponding to the quantum number n = 45. 

FIG. 11. The logarithm of our approximation ( PT| ) (black curve) and the Husimi func- 
tion of the eigenstate n = 30 (gray curve) along the valley of local minima (full line in 
FIG.|^). The curves are projected onto the P-axis. The vertical lines stress the position 
of the Husimi zeroes over the P-axis. The full vertical lines are for zeroes in the principal 
valley of the Husimi function; the dashed vertical lines are for zeroes in the shortest bifur- 
catingva//e?/. The symbols (o) and (o) correspond to the position on the P-axis of the points 
in FIG.| (c) and (d). 

FIG. 12. Idem FICP for the eigenstate n = 45. 
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FIG. 3. 
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